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' 1 INTRODUCTION 



ABSTRACT 

We present a survey of helium abundance indicators derived from a comprehensive 
study of globular cluster photometry in the literature. For each of the three indicators 
used, we have conducted a thorough error analysis, and identified systematic errors 
in the computational procedures. For the population ratio R ~ Nhb /Njigb, we find 
that there is no evidence of a trend with mctallicity, although there appears to be 
real scatter in the values derived. Although this indicator is the one best able to 
provide useful absolute helium abundances, the mean value is 1" ~ 0.20, indicating 
the probable presence of additional systematic error. 

For the magnitude difference from the horizontal branch to the main sequence A 
and the RR Lyrae mass-luminosity exponent A, it is only possible to reliably determine 
relative helium abundances. This is due to continuing uncertainties in the absolute 
metallicity scale for A, and uncertainty in the RR Lyrae temperature scale for A. Both 
indicators imply that the helium abundance is approximately constant as a function 
of [Fc/H]. According to the A indicator, both Oostcrhoff I and II group clusters have 
constant values independent of [Fe/H] and horizontal branch type. In addition, the 
two groups have slopes d log (Pab) /c^pe/H] that are consistent with each other, but 
significantly smaller than the slope for the combined sample. 
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The determination of helium abundances in stars has been 
a long-standing problem in astronomy. For the old stars in 
globular clusters and in the Galactic halo, it is only possi- 
ble to spectroscopically measure the helium abundance in 
a direct way for hot horizontal branch (HB) stars. Studies 
of this kind (e.g. Moehler, Heber, & Durell 1997; Moehler, 
Heber, & Rupprecht 1997) have found abundances signifi- 
cantly higher and lower than the primordial helium abun- 
dance determined from observations of low-metallicity ex- 
tragalactic H II regions (e.g. Yp = 0.234 ±0.002 from Olive, 
Skillman, & Steigman 1997). The interpretation typically 
given for the low helium abundances is that gravitational 
settling has acted preferentially on helium atoms in the at- 
mospheres of these stars, rendering the measurements un- 
usable for determining the initial helium abundance of the 
stars. High helium abundances may be the result of mass 
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loss before the HB phase, which could remove most of the 
hydrogen-rich envelope of the stars. 

Indirect methods must be used to measure the helium 
abundance in earlier phases of stellar evolution. The most 
useful methods are applied to globular clusters because in- 
dividual clusters provide us with large samples of stars with 
seemingly identical compositions and ages. Globular cluster 
stars also seem to be among the oldest that can be found 
in the Galaxy, so that if the helium abundance could be de- 
termined, it would give us a hint of the primordial helium 
abundance. 

There are two other reasons for examining the helium 
abundances of globular clusters. First, it is important to look 
at the overall trend in the helium abundance as a function 
of the cluster metal content. This provides a check of our 
understanding of nucleosynthesis and Galactic evolution. In 
addition, previous studies of RR Lyrae stars have indicated 
that an anticorrelation of helium abundance with metal con- 
tent can explain the Sandage period-shift effect (Sandage 
1982). 

Second, it is important to see if it is possible to de- 
tect helium variations among clusters of similar metallicity, 
since helium remains a plausible candidate for the "second 
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parameter" determining the morphology of the HB. While 
age seems to play a role in some clusters, it appears to be 
inadequate to explain all of the variations seen. In the past, 
clusters have been marked as possibly having high helium 
abundances (e.g., Dickens 1972), but to date there has not 
been definitive evidence of a relative helium abundance dif- 
ference. 

We will examine three different helium abundance indi- 
cators that can be measured from cluster photometry. We do 

this in the hopes that anomalous clusters will appear with 
unusual values in two or more of the parameters, so as to 
provide better evidence for variations. A cursory examina- 
tion of older tabulations of data reveals an additional reason 
for recomputing values for these indicators; the trend of he- 
lium abundance with metallicity appears to depend on the 
method used (Caputo & Castellani 1983). 

The three indicators have been discussed to varying de- 
grees in earlier studies: R (Iben 1968; Buzzoni ct al 1983; 
Caputo, Martinez Roger, & Paez 1987), A (Caputo, Cayrel, 
& Cayrel de Strobel 1983; hereafter, CCC), and A (Carney 
1980; CCC). Photometric studies of globular clusters have 
become rather numerous in the past decade, so that the size 
and quality of the dataset for each parameter can be consid- 
erably increased. Where it has been possible, we have also 
analyzed data for old globular clusters in the Magellanic 
Clouds, and for stars in Local Group dwarf spheroidals with 
old populations. 

In the following sections, wo discuss each parameter in 
turn - their definitions, the datasets, potential errors in mea- 
surement, and calibration. In the final two sections, we com- 
pare helium abundance values derived from the different in- 
dicators, and examine the evidence for trends as a function 
of [Fe/H]. 

2 HELIUM ABUNDANCE INDICATORS 

For all three of the indicators discussed below, higher values 
imply higher helium abundances. 

2.1 The Population Ratio R 

This indicator is simply defined as the ratio R = 
Nhb/Nrgb, where Nhb is the number of horizontal branch 
stars, and Nrgb is the number of RGB stars brighter than 
the luminosity level of the HB. This primarily reflects the 
dependence of the hydrogen-burning shell's progress on the 
hydrogen content of the envelope material being fed into it 
(and to a lessor extent, the change of helium core mass at 
helium flash, which affects the HB luminosity, and hence, 
the HB lifetime). 

Because the ratio is computed using the brightest stars 
in each cluster, it can be calculated for any cluster hav- 
ing deep enough photometry. However, because those stars 
are relatively scarce in globular clusters, it usually requires 
wide-field data, or accurate photometry of the cluster core. 
We have restricted our sample to those clusters with photo- 
metric samples of at least a total of 100 HB and RGB stars. 
In addition, wc require photometry of sufficient quality to 
separate the HB, RGB, and AGB populations readily, and 
to eliminate field stars from the sample when proper motion 
data is not available. 



We have chosen to use the ratio R rather than the simi- 
lar ratio R' = NhbUNrgb + Nagb), which would be more 
easily measured in most cases, because R is more straightfor- 
ward to calculate theoretically. For this reason, the stellar 
sample was sometimes restricted in radius in order to re- 
move portions of the cluster center where image blending 
made identifications of HB, RGB, and HB stars difficult. 

In re-examining the determinations of R, we found that 
the vast majority of earlier calculations used the "average" 
differential bolomctric correction AVbc = Vrgb — Vhb of 
0.15 mag to determine the faint end of the RGB sample. This 
value was used in the original Buzzoni et al. (1983) paper, 
but only because uncertainties of various sorts made the 
correction unimportant. However, the correction is actually 
a function of metallicity, and it is applied to the faint end 
of the RGB sample, which means it can be the source of 
significant systematic error. 

To update AVbc, we have examined the HB models of 
Dorman (1992) in conjunction with the isochrones of Berg- 
busch & VandenBerg (1992; hereafter BV92). These stellar 
models were computed with a consistent set of physics and 
compositions. Although the composition is somewhat out of 
date (it does not include full a-elcmcnt enhancement), the 
differential bolomctric corrections should be satisfactory be- 
cause they are differential in nature, and the oxygen abun- 
dance does not affect the color of the giant branch relative to 
the instability strip. We have determined the corrections as 
a function of [Fe/H], and have derived the following fitting 
formula: 

AVbc = 0.709 + 0.548[M/H] + 0.229[M/H]^ + 0.034[M/H]^ 

Because a-element enhancements influence the position of 
the HB and RGB in the CMD hke a change in [Fe/H] 
(Salaris, Chiefii, & Straniero 1993), they must be taken into 
account when computing [M/H]. To do so, we assumed a 
constant a-clcmcnt enhancement of 0.3 dex for all of the 
globular clusters examined (Carney 1996). [The halo field 
population is believed to show a different abundance pat- 
tern as a function of metallicity (Wheeler, Sneden, & Truran 
1989).] The a-element abundance was taken to contribute 
to the effective metal content of the cluster as 0.7 [a/Fc] 
(Salaris, Chieffi, & Straniero 1993). The contribution does 
not go directly as [a/Fe] because oxygen, the most abundant 
a element, has a relatively high ionization potential. As a re- 
sult, it does not contribute as significantly to the opacity in 
the envelope of RGB stars as it does for higher temperature 
HB stars. 

Error in observed R values comes from two sources: er- 
ror in the numbers of HB and RGB stars due to misidenti- 
fication or Poissori fluctuations, and error in the determina- 
tion of the faint limit for the RGB sample (resulting from 
the determination of the magnitude level of the HB, or the 
metallicity uncertainty that affects the size of the differential 
bolomctric correction) . 

The magnitude level of the HB was taken to be the 
average magnitude of stars in the instability strip. For clus- 
ters with many well-studied RR Lyrae stars, this could be 
taken from the average magnitude of the variables. Often 
though the RR Lyrae stars were not observed frequently 
enough to derive good average magnitudes. If the cluster 
had populated red and blue edges of the intability strip, a 
linear interpolation between the edges was made. In cases 
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where the cluster had only a populated red or blue HB, a 
theoretical correction was made using the models of Dor- 
man (1992) to logj^Q Teff = 3.85 using the cluster reddening 
and a well-determined point on the populated portion of 
the HB. Dorman's scaled-solar abundance HB models were 
used with a correction for enhanced a-element abundances 
as described above. 

An examination of the table indicates that the quoted 
errors for the HB magnitudes for the clusters are in general 
less than 0.05 mag. This is because we are not interested in 
the absolute value of Vhb for any of the clusters - only the 
value for the sample used. It is clearly unnecessary for the 
photometric calibration of the sample to be perfect, since 
we only require that the relative photometry of the HB and 
RGB be good enough that the faint limit of the RGB sample 
can be determined well. For this reason, the tabulated values 
of Vhb should not be taken as being good absolute values. 
This column is provided to make verification of the results 
easier. 

The error due to the differential bolometric correction 
is illustrative, so we present a semi-analytic model for its 
contribution to the error. The important quantity to com- 
pute is ONrgb /OVrgb ■ This can be derived assuming that 
i) the differential luminosity function (LF) of the RGB has 
a constant slope a = d(logj^Q N)/dV, and ii) the number of 
red giant branch stars goes to zero at the tip of the RGB 
(TRGB). So, we use logio(iV) = a{V - Vtrgb), and find 



N, 



-(10- 



'(.VRGB—yTRGB) 



In 10 

where No is the normalization of the luminosity function. 
From this we find 

^Qa{VfiQB —Vtrgb) 



dNuGB 

OVrgb 



= a In 10 



IQo'iVRGB -Vtrgb) 



Nrgb- 



The fraction in brackets is close to one, and from Bergbusch 
& VandenBerg (1992) LFs we find a = 0.33 nearly indepen- 
dent of composition. From this we find that for each tenth 
of a magnitude added to the faint limit of the RGB sample, 
the number of RGB stars increases by about 9%. (Alter- 
nately, making the faint limit fainter by one magnitude for 
a given cluster area increases the total number of red giants 
by approximately 85%.) This demonstrates the importance 
of accurate values for the difi'erential bolometric correction 
AVbc- 

The above derivative affects the error in Nrgb through 
errors in metallicity and in the determination of the HB 
magnitude. Because of the metallicity dependence of the 
differential bolometric correction, helium determinations for 
high metallicity clusters will be relatively more uncertain. A 
polynomial fit to the derivative gives 

dln NRGB ^ Q g^gg ^ 0.1084[M/H1 + 0.0257[M/H1^ 
oVrgb 

We note that the error from this source goes as Nrgb, 

1/2 

whereas the Poisson counting errors go as N^q^. 

We have used [Fe/H] values from the tabulation of Djor- 
govski (1993) in computing the values of AVsc for the clus- 
ters. Djorgovski's values fairly closely follow the values of 
Zinn & West (1984; hereafter, ZW). We have also consid- 
ered an alternative metallicity scale since this affects the 
differential bolometric corrections. Using metallicities from 
Carretta & Gratton (1996), we recomputed the R values 



Figure 1. A comparison of _R values derived assuming metallic- 
ities from the scale of Zinn & West (1984; ZW) and Carretta & 
Gratton (1996; CG). 



for 16 clusters from their sample. The results are shown in 
Fig. |l]. For the most part, the values were little affected by 
the changes in the scale. 

Computed R values are insensitive to factors that 
merely shift the RGB in temperature or color (age, for exam- 
ple) . Factors that cause changes in the absolute brightness of 
the HB (such as mean mass of the helium cores of stars, and 
CNO abundances) will cause systematic errors in measured 
values for individual clusters. 

Our sample is composed of 42 Galactic globular clus- 
ters (GGCs) and 5 globulars in the Large Magellanic Cloud 
(LMC), and is presented in Table ^. The listings are grouped 
in bins of approximately 0.2 in [Fe/H]. References carrying 
a "(PM)" designation were used to remove field stars from 
the sample using proper motion information. There tend to 
be few clusters at the very metal-rich end due to a combi- 
nation of substantial field star contamination and confusion 
between RGB and red HB cluster stars. 

The variation of the R values as a function of HB type 
is shown in Fig. |^. There may be a slight decrease in the 
average R for clusters with bluer HB type. Several clusters 
with the bluest morphologies have unusually high R values: 
NGC 6752, NGC 7099 (M30), NGC 6341 (M92), NGC 6218 
(M12), and NGC 6681 (M70). However, this is not universal 
— there are more clusters in the same range of Rhb having 
R values that are closer to the Galactic average. There is 
not a clear reason for this difference: well-measured clusters 
are found in both groups, and a metallicity difference does 
not appear to be present. 

An important effect of the inclusion of the metallicity 
dependence of AVbc is the improved agreement between 
the R values of several of the most metal-rich clusters and 
the metal-poor clusters. In a similar fashion, the increased 
AVbc values will also significantly reduce the estimates of 
the helium abundances of Galactic bulge fields {Y — 0.28 ± 
0.02 according to Minniti 1995) because of the high mean 
metallicity of those stars. 

2.2 The MS-HB Magnitude DifFerence A 

The indicator A (CCC) is defined simply as the magnitude 
difference between the MS at {B — V)o = 0.7 and the HB at 
the instability strip. CCC originally defined the HB point to 
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Table 1. Data for the bright-star population ratio R 



ID (NGC/IC) 




Vhb 


AVbc Nrgb 


J^AGB 


R 


References 
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10 
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23 
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7099 (M30) 


202 
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0.27 
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Sandquist ct al 1998, 



Yanny et al 1994 
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Table 1 - continued 



ID (NGC/IC) 


Nhb 


Vhb 




I^RGB 


Nagb 


R 


References 








Magellanic Clouds Clusters 






1466 


91 


19.33 ±0.02 


0.28 


96 


36= 


0.95 ±0.14 


Walker 1992c 


1841 


165 


19.35 ±0.03 


0.26 


158 


32 


1.04 ±0.12 


Brocato et al 1996c 


2210 


67 


19.10 ±0.03" 


0.27 


70 


15 


0.96 ±0.17 


Brocato et al 1996c 


2257 


164 


19.02 ±0.03 


0.28 


138 


24 


1.19 ±0.14 


Testa et al 1995 


Reticulum 


46 


19.07 ±0.01 


0.29 


54 


12 


0.85 ±0.17 


Walker 1992b 



Notes: " A theoretical correction has been applied to compute the level of the instability strip. 
^ The possibility of incompleteness exists at the faint end of the sample. 
The number may be affected field star contamination. 



Figure 2. Globular cluster R values versus their horizontal 
branch types (Rhb = {B — R) / (B + V + R)). • indicates clusters 
having greater than 200 observed HB stars, o indicates clusters 
having between 100 and 200, * indicates clusters with less than 
50, and A indicates LMC clusters. 



be at the blue edge of the instability strip, but we have cho- 
sen to revise this definition to make it more easily calculable 
theoretically and observationally. This also reduces the color 
difference between the HB and MS points, reducing possible 
systematic effects from photometric calibration. The various 
sensitivities of the indicator are not significantly changed by 
the revision to the definition. 

Increases in envelope helium abundance infiuence this 
indicator by increasing the luminosity of the HB (via the 
strength of hydrogen shell-burning during that phase) and 
by increasing the effective temperature and the luminosity of 
the MS through a decrease in the envelope opacity. The net 
effect of the luminosity and temperature changes on the MS 
is to make it fainter at a given color. As an indicator, A has 
the advantage of being sensitive to the helium abundance 
{dA/dY = 5.8 mag; CCC), and the disadvantages of hav- 
ing a definite metallicity dependence (9A/9[Fe/H] « —0.5 
mag / dex) and of requiring photometry from the HB to 
well below the MS turnoff. There is an additional disadvan- 
tage in requiring the knowledge of the cluster's reddening 
to determine Vo.r- However, A has no dependence on age, 
since the chosen MS point reaches unevolved stars. We can 
calculate theoretical values of A using a self-consistent set 
of isochrones (BV92) and HB models (Dorman 1992) hav- 
ing Y « 0.236. We have fitted the following polynomial to 
the data after making small corrections for different initial 
helium content in the theoretical models: 



A = 4.268 - 2.1295[M/H] - 0.7938([M/H])^ 

-0.1173([M/H])^ 

As can be seen from Table ^, there are only 20 clusters 
with enough information to compute a value of A. We in- 
cluded clusters for which Vhb and V0.7 were derived from 
different studies, in spite of the possibility of zero point dif- 
ferences. The error budget for A was computed using: 

a^iA)^a\VHB)+a\Vo.7) 

+ (^(f^)V(i5(B-l/)) + a^(0) 

where the last two terms account for reddening and zero 
point uncertainties respectively. As was discussed in § 2.1, 
we do not need to worry about how well calibrated the data 
are in an absolute sense, as long as there are not significant 
systematic errors in the relative photometry for the MS and 
HB (for instance, a nonlinearity in one of the CCDs). 

In a global sense, the degree of agreement between the 
observed values and theoretical expectations is primarily a 
function of the metallicity scale chosen. In Fig. |3[ we show 
the A values using metallicities from three different studies: 
ZW, Carretta & Gratton (1996; hereafter CG), and Rut- 
ledge, Hesser, & Stetson (1997; hereafter RSD). The ZW 
comparison agrees well with the theory, although the scatter 
in values at a given metallicity is large. The high-resolution 
spectroscopy study of CG reduces the scatter considerably, 
although there seems to be a constant offset between the 
observed and theoretical values. RSD calibrated their Ca 
II triplet measurements to both the ZW and CG scales. In 
Fig. H, we plot the values using RSD's CG scale. Again, 
the scatter is fairly low, but the shape of the curve traced 
by the observed values is very different from theoretical ex- 
pectations. The situation is almost identical for their ZW- 
calibrated values. 

How are we to judge the merits of the different compar- 
isons? First we must keep in mind that the absolute metal- 
licity scale is uncertain at the 0.2 dex level (RSD), as has 
been assumed throughout the paper. The relative rankings 
are to be trusted more {a ~ 0.1 dex), making the shape 
of the A comparison most reliable. So, horizontal shifts of 
the curve by 0.2 dex should not be considered unreasonable. 
This means that the comparisons using both the ZW and CG 
metallicity data are consistent with the theoretical values. 
The RSD scales result in large discrepancies, particularly at 
the metal-rich end. The fact that the shape is different in 
both their ZW and CG calibrations implies that it probably 
results from details or assumptions of their technique. RSD 
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Table 2. Data for the magnitude difference A 

ID (NGC) Vhb E(B - V) Vo.7 A HB Reference MS Reference 



6838 {M71) 


14, 


.41 


± 


0. 


,03" 


0, 


.28 ± 


0. 


,02 


19, 


.80 


± 


0, 


.12 


5, 


.39 


± 


0, 


.12 


Hodder et al 1992 


Hodder et al 1992 


104 (47 Tuc) 


14, 


.02 


± 


0. 


,05" 


0, 


.04 ± 


0. 


,01 


19, 


.43 


± 


0, 


.09 


5, 


.41 


± 


0, 


.07 


Hesser et al 1987 


Hesser et al 1987 


6171 


(M1U7) 


15, 


.62 


± 


0. 


,02 


0, 


.31 ± 


0. 


,02 


20, 


.81 


± 


0, 


.13 


5, 


.19 


± 


0, 


.13 


Ferraro et al 1991 


Ferraro et al 1991 


362 




15, 


.42 


± 


0. 


,02" 


0, 


.04 ± 


0. 


,01 


21, 


.36 


± 


0, 


.07 


5, 


.94 


± 


0, 


.09 


Harris 1982 


Bolte 1987 


1261 




16, 


.72 


± 


0. 


,03 


0, 


.00 ± 


0. 


,02 


22, 


.32 


± 


0, 


.11 


5, 


.60 


± 


0, 


.12 


Alcaino et al 1992 


Alcaino et al 1992 


1851 




16, 


.15 


± 


0. 


,02 


0, 


.02 ± 


0. 


,02 


21, 


.60 


± 


0, 


.09 


5, 


.45 


± 


0, 


.09 


Walker 1992a 


Walker 1992a 


6121 


(M4) 


13, 


.50 




0. 


.09 


0, 


.37 ± 


0. 


.01 


18. 


.98 


± 


0. 


.06 


5, 


.48 




0, 


.11 


Kanatas et al 1995 


Richer & Fahlman 1984 














































Alcaino et al 1988 


288 




15, 


.26 




0. 


,04" 


0, 


.03 ± 


0. 


,01 


21. 


.17 




0, 


.07 


5, 


.91 




0, 


.10 


Bolte 1992 


Bolte 1992 


5139 


(u) Cen) 


14, 


.56 


± 


0. 


,01 


0, 


.15 ± 


0. 


,02 


20, 


.60 


± 


0, 


.13 


6, 


.04 


± 


0, 


.16 


Kaluzny et al 1997 


Alcaino & Liller 1987 


5904 


(M5) 


15, 


.09 


± 


0. 


,02 


0, 


.03 ± 


0. 


,01 


20, 


.88 


± 


0, 


.06 


5, 


.78 


± 


0, 


.06 


Sandquist et al 1996 


Sandquist et al 1996 


6752 




13, 


.70 


± 


0. 


,05" 


0, 


.04 ± 


0. 


,01 


19, 


.69 


± 


0, 


.07 


5, 


.99 


± 


0, 


.09 


Penny & Dickens 1986 


Penny & Dickens 1986 


5272 


{M3) 


15 


.66 


± 


0. 


03 


0, 


.00 ± 


0. 


.01 


21, 


.50 


± 


0, 


.05 


5 


.84 




0, 


.06 


Ferraro et al 1997 


Ferraro et al 1997 


6218 


{M12) 


14, 


.90 


± 


0. 


,05 


0, 


.22 ± 


0. 


,04 


20, 


.43 


± 


0, 


.28 


5, 


.53 


± 


0, 


.30 


Racine 1971 


Sato et al 1989 


6254 


(MID) 


14, 


.62 


± 


0. 


,05" 


0, 


.27 ± 


0. 


,03 


20, 


.68 


± 


0, 


.12 


6, 


.06 


± 


0, 


.13 


Hurley et al 1989 


Hurley et al 1989 


6397 




12, 


.97 


± 


0. 


,05" 


0, 


.18 ± 


0. 


,02 


19, 


.27 


± 


0, 


.10 


6, 


.30 


± 


0, 


.11 


Alcaino et al 1987 


Alcaino et al 1987 


4590 


{M68) 


15, 


.64 


± 


0. 


,01 


0, 


.07 ± 


0. 


,01 


22, 


.00 


± 


0, 


.05 


6, 


.36 


± 


0, 


.05 


Walker 1994 


McClure et al 1987 


6341 


{M92) 


15, 


.10 


± 


0. 


,01 


0, 


.02 ± 


0. 


,01 


21, 


.40 


± 


0, 


.05 


6, 


.30 


± 


0, 


.05 


Carney et al. 1992b 


Stetson & Harris 1988 


7078 


{M15) 


15, 


.81 


± 


0. 


,01 


0, 


.10 ± 


0. 


,01 


22, 


.15 


± 


0, 


.05 


6, 


.34 


± 


0, 


.05 


Silbermann & Smith 1995b 


Durrell & Harris 1993 


7099 


{M30) 


15, 


.10 


± 


0. 


,02 


0, 


.05 ± 


0. 


,02 


21, 


.61 


± 


0, 


.10 


6, 


.52 


± 


0, 


.10 


Bergbusch 1996 


Bergbusch 1996 



Note: " Vhb has been corrected to the instability strip. 



ods) as it is for the CG and RSD studies since they both 
applied a single method in a uniform way. We have chosen 
to use the CG scale because of the general agreement of 
the observational and theoretical curve shapes. The helium 
abundance for each of the clusters was computed relative to 
the theoretical values given the derivative dA/dY (CCC). 
In Figs. ^ and ^ we plot the calculated 5Y values as a 
function of metallicity and HB morphology, respectively. It 
is unsettling that the clusters with |-R//fl| < 0.8 (and well- 
defined blue and red edges to the instability strip) tend to 
have higher SY values. NGC 288 and M30 stand out slightly 
as having high values among the clusters with the bluest 
morphologies, and NGC 362, M71, and 47 Tuc stand out as 
high for clusters with red morphologies. This may indicate 
that the level of the horizontal branch is not being mea- 
sured properly. However, the magnitude error would have 
to be between about 0.1 and 0.25 mag to bring them back 
into agreement, which seems overly large. If the reddest HB 
clusters can be explained via errors in their metallicities, 
there might be slight evidence for increased helium abun- 
dance with bluer HB type. We will discuss discrepant indi- 
vidual clusters in § |^. 

2.3 The RR Lyrae Mass-Luminosity Exponent A 

The indicator A (CCC) is related to the mass-luminosity 
relationship for stars inside the instability strip: 

A = log(L/i0) - 0.81 log(M/M0). 

A is dependent on helium because increased helium can 
both increase the luminosity of an RR Lyrae, as well as the 



Figure 3. The A indicator as as function of [Fe/H]. The three 
panels use metallicities taken from the studies cited. The Rut- 
ledge et al. panel uses their metallicities calibrated to Carretta &; 
Gratton's scale. 

discuss the issue of their Ca measurements as indicators of 
[Fe/H] at length, and we refer readers to their § 6. 

Because the absolute metallicity scale remains uncer- 
tain at about the 0.2 dex level, it is currently impossible to 
determine low-error absolute helium abundances using A. 
At the metal-rich end in particular, the absolute metallic- 
ity uncertainty causes large uncertainties in absolute helium 
abundance. The relative rankings are to be trusted more if 
the relative metal abundances for the globular clusters are 
good. This is not as much the case for the ZW scale (since it 
was a weighted average of determinations by different meth- 
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Figure 4. SY{A) values as a function of Rhb (Lee, Demarque, 
& Zinn 1994). 

mean mass of stars occupying the instability strip. While 
A has a relatively small sensitivity to helium abundance 
{dA/dY = 1.4; Sandage 1990b, Bencivenni et al. 1991), star 
tistical errors are generally small for clusters with a fair num- 
ber of RR Lyracs (a a ~ 0.01). Potentially A could provide 
helium abundances with the best precision (cry ^ 0.007) of 
the three indicators we have considered. 

{A) can be computed using the period relation of van 
Albada & Baker (1971): 

log P = 11.497 + 0.84^ - 3.481 log Te//. 

One of the difficulties in using this relation is the computa- 
tion of realistic effective temperatures for the cluster vari- 
ables from readily observable quantities. As will be shown 
below, the uncertainty in the absolute temperature values 
results primarily in systematic shifts in {A) values. As a 
result, (A) can only be realistically considered a relative in- 
dicator of helium abundance at this time. 

One source of uncertainty in the temperature scale re- 
lates to the calibration samples of field RR Lyraes. Aside 
from uncertainty in the model atmospheres used to cali- 
brate the colors, the optimal choice of color continues to 
be debated. {V — K) has been recommended by many au- 
thors (Liu & Janes 1990; Jones et al. 1992; Carney, Storm, 
& Jones 1992a; Fernley 1993) over [B — V) due to evidence 
of shock- wave effects on B magnitudes near maximum light. 
However, McNamara (1997) showed that the temperature 
calibration from {V — K) disagreed systematically as a func- 
tion of period with those from {B — V),{h — y),{y — R), and 
{V — I) when 1994 Kurucz model atmospheres were used. 

In recent years, studios of globular cluster RR Lyraes 
have turned to the use of quantities like the period log P and 
blue amplitude As to derive temperatures so as to avoid 
systematic errors resulting from reddenings. However, there 
is still uncertainty in the temperature zero-point from the 
model atmospheres used to calibrate the temperature of field 
RR Lyraes, and from the differences in temperatures derived 
from colors using different filter combinations. 

We have chosen to recalibrate the temperature relation 
of Catelan (1998) for RR Lyrae stars of type ab (with Ab 
and [Fe/H] the only variables) to temperatures derived us- 
ing (B — V) colors. Catelan, Sweigart, & Borissova (1998a) 
point out that the use of a log P term in the determination 
of effective temperatures (as in equation 16 of Carney et al. 



(1992a)) tends to cause luminosity differences among RR 
Lyraes to be translated into temperature differences, erro- 
neously reducing the scatter in the P — T^q plane. 

Our decision to use {B — V) color temperatures for the 
calibrating sample of field RR Lyraes is based on the find- 
ings of McNamara (1997), and Kovacs & Jurcsik (1997). 
As mentioned above, McNamaxa found that temperatures 
derived from {V — K) deviated systematically from several 
other commonly-used colors as a function of period. It has 
been known for a long time that Mk correlates strongly 
with logP (e.g. Longmore et al. 1990). Because the K band 
is on the Ray leigh- Jeans portion of the blackbody curve for 
RR Lyraes, it is not very sensitive to temperature [B and 
V fall near the maximum). Kovacs & Jurcsik's examination 
of Mv : Ml,, , and Mk as a function of Fourier light curve 
parameters for globular cluster variables indicates that the 
period dependence is important for each of the filters, but is 
larger for redder filters. As a result, filter combinations with 
longer wavelength baselines have larger period dependences. 
In particular, {V — K) is predicted to have a period de- 
pendence that can explain the trend McNamara sees, while 
(V—Ic) has a dependence that is over twice as small, and the 
dependence for {B — V) is over ten times smaller. Because 
period depends significantly on the stellar luminosity in ad- 
dition to temperature (see van Albada & Baker's pulsation 
equation), a large period dependence in the color is likely to 
be a systematic problem, in agreement with the Eissertion of 
Catelan et al. (1998a). 

So, using temperatures from McNamara (1997), pulsa- 
tional amplitudes from Blanco (1992), and metallicities from 
Layden et al. (1996) for field RR Lyraes, we found the relar 
tion: 

e^g = (0.776 ± 0.012) - (0.030 ± 0.008)As 

-(0.008 ±0.003) [Fe/H]. 

The fit has a multiple correlation coefficient r = 0.925, and 
an rms deviation of 40 K from the fit. In addition, the resid- 
uals show no correlation with log P. To check the effect 
that different compositions for field and cluster RR Lyraes 
could have, we redid the fit using only variables with [Fe/H] 
< —1.0. The fit (using 15 stars) was 

Qeq = (0.786 ± 0.020) - (0.039 ± 0.011)Ab 

-(0.008 ±0.009) [Fe/H]. 

We have chosen to use this second calculation for the calcu- 
lations in this paper, though it makes only small changes to 
the (A) values. 

In computing {A} values for cluster variables, we have 
chosen not to include variables which are known to exhibit 
the Blashko effect, as this causes changes in pulsation am- 
plitude (and in our computed temperatures) from cycle to 
cycle. Those variables were included in calculations of the 
average period though, since the period of pulsation is not 
affected. Error analysis was carried out on A values for each 
variable. However, we find that the error in individual A 
values is not very significant compared with the scatter in A 
values for variables within a cluster. From an examination 
of histograms of A values for the most populated clusters, 
the distributions appear to be Gaussian. 

We have found 50 Milky Way globular clusters having 
good data on at least 2 RRab variables, for a total of 974 
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Figure 5. A comparison of (A) values from the present study 
witli (A(12)) values from Caputo & De Santis (1992). 

stars. We have also analyzed 8 old Magellanic Cloud clus- 
ters (108 variables), and 4 Local Group dwarf spheroidal 
galaxies (214 variables), as shown in Table ^ Only 25 of the 
Milky Way clusters (and 5 in the Magellanic Clouds) have 

10 or more RRab stars with good data. The table includes 
only those clusters for which there is at least one RRab star 
with period and amplitude measurements. For the clusters 
examined, we give the mean period of the RRab stars in 
column (3), along with the error in the mean. The average 
mass-luminosity exponent {A) is computed from the average 
of the values for the stars that have Ab values. When this 
is different from the number that have periods, the number 
is given in parentheses in column (2) . Fig. ^ plots our data 
against those calculated according to the method of Caputo 
& De Santis (1992), which employed a different (but also 
reddening-independent) method for computing A. The pri- 
mary difference between the studies is in zero-points, which 
is probably due to the differences in the temperature cali- 
brations. There may also be a slope difference at high values 
of (A) . 

Caputo & De Santis found that their (A) values showed 
sensitivity to the HB type — clusters with blue HBs {{B — 
R)/{B + V + R) > 0.7) tended to have significantly higher 
values for {A) . Synthetic HB models predict that only very 
evolved HB stars are found in the instability strip for clus- 
ters with very blue HB populations. In Fig. we plot our 
{A) values versus horizontal branch type Rhb- As Caputo 
& De Santis found, there is little scatter among Ooster- 
hoff group I clusters with a few exceptions (most notably 
NGC 6229 on the low end). None of the Oosterhoff group 

11 clusters have {A) values consistent with the overwhelm- 
ing majority of Oo I clusters. (Although Ruprecht 106 has 
a {Pab) closer to those of the Oosterhoff II group, its lack 
of RRc variables and its low (A) value indicate that it could 
be an Oo I cluster.) M14 and M28, the two Oo I clusters 
with the bluest HB morphologies, have {A) values that place 
them at the low end of the range populated by OosterhofT 
II clusters. 

There are a number of systems that have (A) values 
that are consistent with being in the Oosterhoff II group 
and HB types that are redder than the majority of Oo II 
clusters. This is partly an incarnation of the second param- 
eter effect in HB morphology. However, it seems to rule out 
the possibility that Oo II RR Lyraes are simply the result of 



Figure 6. The variation (A) writh horizontal branch type {Rhb = 
{B — R)/{B + V + R)). Solid symbols indicate clusters with 7 or 
more RRab variable stars. (All of the dwarf spheroidals tabulated 
have samples this large.) 

the instability strip being populated solely by stars that are 
evolving toward the AGB. The instability strip is too well 
populated in systems like the Ursa Minor dSph and Galactic 
globular clusters like M68 and NGC 5466, all of which have 
populations of red HB stars in addition to the variables. 

For this reason, we have re-examined the dependence 
of (A) with [Fe/H] in Fig. |. Many previous studies (e.g. 
Sandage 1990b, Caputo & De Santis 1992) have noted the 
apparent linear correlation within the total globular cluster 
sample. If this is truly the case, then the two Oosterhoff 
groups should also show the same linear relation separately. 
We derive the following relations; 

A = (-0.032 ± 0.008)[Fe/H] + (1.821 ± 0.012) 

for the 21 Oo I clusters having 7 or more RRab stars (ex- 
cluding NGC 6229 as an extreme outlier) having —1.69 < 
[Fe/H] < -0.99, and 

A = (0.001 ± 0.018) [Fe/H] + (1.929 ± 0.036) 

for 12 Oo II clusters with 7 or more RRab stars (excluding cu 
Cen due to its metallicity spread) having —2.22 < [Fe/H] < 
— 1.58. If Ruprecht 106 is removed from the Oo I sample (it 
appears to be a young cluster, and it has an unusually low 
dispersion in the A values of its variables), we find 

A = (0.003 ± 0.010)[Fe/H] + (1.868 ± 0.015). 

The two groups individually have slopes that are signifi- 
cantly shallower than derived from the union of the two 
samples (— 0.088±0.006), and both are consistent with zero. 
Both slopes are also significantly smaller than predictions 
from synthetic HB computations (—0.027; Lee, Demarque, 
& Zinn 1990). For the Oo I clusters the small slope is not 
surprising since the {A) values are not predicted to have a 
dependence on HB type for Rhb < 0.7. The Oo I clusters 
M14 and M28 are the only Galactic globular clusters (with 
40 and 9 RRab variables respectively) that can be said to fall 
in the gap between the two groups. Both have very blue HB 
morphologies, indicating that the relatively high tempera- 
tures of the variables are affecting the {A) values (Caputo et 
al. 1993). More unexpected is the essentially constant value 
found for the Oo II clusters which are expected to be much 
more sensitive to HB type (Caputo & De Santis 1992). 
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Table 3. Data for the RRab mass-luminosity relation (A) 



ID (NGC/IC) 


Nab 


(Pab) 








(A) 


References 














Milky Way Clusters 


in/i /a^ '~r\tr^\ 
iU4 (4( lucj 


1 


0.737 








2.039 


Carney et al 1993 


6388 


4 


0.746 


-|- 


0. 


.071 


2.039 ± 0.044 


Hazen & Hesser 1986, 


6171 (M107) 


15 


0.537 


± 


0, 


.018 


1.853 ±0.013 


Dickens 1970 


6366 


1 


0.508 








1.852 


Clement 1997 


6362 


15 


0.556 


± 


0, 


.018 


1.874 ±0.012 


Clement et al 1995a 


6638 


1 


0.666 








1.931 


Clement 1997 


6712 


7 


0.557 


± 


0, 


.022 


1.875 ±0.012 


Sandage et al 1966 


6723 


23(20) 


0.536 


± 


0, 


.014 


1.862 ±0.012 


Menzies 1974 



362 


7 


0, 


,542 


± 


0, 


,023 


1261 


13 


0. 


,555 


± 


0, 


,013 


1851 


21(18) 


0. 


,571 


± 


0, 


.014 


2808 


1 


0. 


,539 








6121 (M4) 


30(25) 


0. 


,538 


± 


0, 


,010 


6266 {M62) 


62(40) 


0, 


,544 


± 


0, 


,009 


6402 (M14) 


40 


0, 


,572 


± 


0, 


,012 


6717 (Pal 9) 


1 


0, 


,575 








6864 (M75) 


3 


0, 


,531 


± 


0, 


,030 



1 


.860 


± 


0, 


.017 


Sawyer Hogg 1973 


1, 


.870 


± 


0, 


.012 


Wehlau & Demers 1977 


1 


.873 


± 


0, 


.014 


Walker 1998 


1 


.857 








Clement & Hazen 1989 


1 


.849 


± 





.009 


Sujarkova & Shugarov 1981 


1 


.851 




0, 


,009 


Sawyer Hogg 1973 


1 


.886 


± 


0, 


,009 


Wehlau & Froclich 1994 


1 


.895 








Clement 1997 


1 


.842 


± 


0, 


,022 


Pinto et al 1982 



288 


1 


0. 


,679 








1, 


,972 








4499 


63(59) 


0. 


,580 


± 


0, 


.008 


1, 


.879 


± 


0. 


,006 


5139 {uj Cen) 


83(74) 


0. 


,646 


± 


0. 


.011 


1, 


.929 


± 


0. 


,009 


5904 (M5) 


84(80) 


0. 


,552 


± 


0, 


,008 


1, 


,863 


± 


0. 


,006 


6229 


11 


0, 


,527 


± 


0, 


,013 


1, 


,818 


± 


0, 


Oil 


6235 


2 


0, 


,600 


± 


0, 


,016 


1, 


,927 


± 


0, 


,010 


6284 


6 


0. 


,588 


± 


0, 


,023 


1, 


,909 


± 


0, 


,012 


6544 


1 


0, 


,57 








1, 


,90 








6558 


6 


0, 


,556 


± 


0, 


,044 


1, 


,899 


± 


0, 


,032 


6584 


34 


0. 


,560 


± 


0, 


Oil 


1, 


,876 


± 


0. 


,007 


6626 (M28) 


9 


0. 


,586 


± 


0, 


,022 


1, 


.895 


± 


0. 


,020 


6681 (M70) 


1 


0. 


,564 








1, 


.911 








6715 (M54) 


27 


0. 


,553 


± 


0, 


,012 


1, 


,856 


± 


0. 


Oil 


6934 


45 


0, 


,561 


± 


0, 


,012 


1, 


,853 


± 


0, 


,008 


6981 (M72) 


25 


0, 


,559 


± 


0, 


,009 


1, 


,851 


± 


0, 


,008 


7006 


53(50) 


0, 


,566 


± 


0, 


,005 


1, 


,861 


± 


0, 


,005 


7089 (M2) 


13 


0, 


,623 


± 


0, 


,020 


1, 


,925 


± 


0, 


,017 


7492 


1 


0. 


,805 








2, 


.015 
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Sawyer Hogg 1973 

Clement 1997 

Clement 1997 

Hazen 1993 

Hazen 1996 

Millis & Liller 1980 

Wehlau & Butterworth 1990 

Liller 1983 

Sawyer Hogg 1973 

Sawyer Hogg & Wehlau 1980 

Dickens & Flinn 1972 

Wehlau et al 1992 

Sawyer Hogg 1973 

Clement 1997 



1904 


2 


0. 


,685 


± 


0. 


.051 


1. 


.952 


± 


3201 


72(66) 


0. 


,555 


± 


0, 


.005 


1, 


.860 


± 


5272 (M3) 


148(92) 


0. 


,553 


± 


0, 


,006 


1, 


,865 


± 


5286 


10 


0, 


,614 


± 


0, 


,019 


1, 


,915 


± 


5897 


4 


0, 


,735 


± 


0, 


,094 


1, 


,956 


± 


5986 


7 


0, 


,650 


± 


0, 


,035 


1, 


,945 


± 


6093 (M80) 


4 


0, 


,651 


± 


0, 


,012 


1, 


,931 


± 


6139 


3 


0, 


,537 


± 


0, 


,113 


1, 


,893 


± 


6273 (M19) 


1 


0. 


,507 








1, 


,838 




6333 (M9) 


8 


0. 


,622 


± 


0. 


.018 


1. 


.940 


± 


6656 (M22) 


9 


0. 


,640 


± 


0, 


.015 


1, 


.922 


± 


Palomar 13 


4 


0. 


,572 


± 


0, 


.012 


1, 


.879 


± 


Rup 106 


12(10) 


0. 


,616 


± 


0, 


.006 


1, 


.883 


± 



0.027 Sawyer Hogg 1973 

0.004 Cacciari 1984 

0.004 Sandage 1959 

0.016 Liller & Lichten 1978a 

0.086 Wehlau 1990 

0.019 Liller & Lichten 1978b 

0.013 Wehlau et al 1990 

0.059 Hazen 1991 

Clement & Sawyer Hogg 1978 

0.010 Clement et al 1984, Clement & Shelton 1996 

0.014 Wehlau & Sawyer Hogg 1978 

0.006 Sawyer Hogg 1973 

0.002 Kaluzny et al 1995 



2298 1 

4147 4 

4833 7 

5634 3 

5824 7 

6293 2 



0, 


,640 








1 


.924 








0, 


,531 


± 


0, 


,029 


1 


.827 


± 


0, 


,031 


0, 


,685 


± 


0, 


,031 


1 


.951 


± 


0, 


,020 


0. 


,621 


± 


0, 


,019 


1, 


.925 


± 


0. 


,012 


0. 


,624 


± 


0. 


.008 


1 


.920 


± 


0. 


,006 


0. 


,600 


± 


0, 


.006 


1 


.933 


± 


0. 


,008 



Clement et al 1995b 
Ncwburn 1957 
Demers & Wehlau 1977 
Liller & Sawyer Hogg 1976 
Clement 1997 
Clement et al 1982 
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Table 3 - continued 



ID (NGC/IC) 


Nab 


x-* 0,0/ 




(A) 


References 


2419 


25 


0.656 ±0, 


.007 


1.933 ±0.006 


Pinto & Rosino 1977 


4590 (M68) 


11(9) 


0.627 ±0, 


.016 


1.922 ±0.014 


Walker 1994 


5024 {M53) 


24(23) 


0.646 ± 0, 


.014 


1.929 ±0.007 


Goranskij 1976 


5053 


5 


0.672 ± 0, 


.026 


1.906 ±0.020 


Nemec et al 1995 


5466 


11(9) 


0.637 ±0, 


.023 


1.921 ±0.018 


Sawyer Hogg 1973 


6341 (M92) 


9(6) 


0.628 ± 0, 


.017 


1.907 ±0.022 


Carney et al 1992b 


6426 


4 


0.665 ± 0, 


.020 


1.926 ±0.015 


Sawyer Hogg 1973 


7078 (MIS) 


31(25) 


0.640 ± 0, 


.010 


1.918 ±0.009 


Silbermann & Smith 1995, Sandagc 1990a 


7099 (MSO) 


3 


0.698 ± 0, 


.026 


1.965 ±0.029 


Sawyer Hogg 1973 








Magellanic Clouds Clusters 


121 


4 


0.548 ± 0, 


.031 


1.845 ±0.026 


Walker & Mack 1988a 


1466 


25(21) 


0.581 ±0, 


.013 


1.876 ±0.008 


Walker 1992c 


1786 


5 


0.677 ±0, 


.042 


1.942 ±0.022 


Walker & Mack 1988b 


1835 


21(18) 


0.587 ±0, 


.017 


1.899 ±0.015 


Walker 1993 


1841 


17(15) 


0.676 ± 0, 


.013 


1.934 ± 0.009 


Walker 1990 


2210 


14(9) 


0.585 ± 0, 


.018 


1.879 ±0.020 


Hazen & Nemec 1992, Reid & Freedman 1! 


2257 


18(15) 


0.592 ±0, 


.018 


1.886 ±0.014 


Nemec et al 1985, Walker 1989 


Reticulum 


22(21) 


0.552 ±0, 


.012 


1.855 ±0.007 


Walker 1992b 










Dwarf Spheroidals 


Carina 


49 


0.620 ±0, 


.006 


1.930 ±0.005 


Saha et al 1986 


Draco 


115(98) 


0.613 ±0, 


.004 


1.888 ±0.003 


Nemec 1985 


Sextans 


26(20) 


0.606 ± 0, 


.010 


1.894 ± 0.009 


Mateo et al 1995 


Ursa Minor 


47 


0.636 ±0, 


.009 


1.922 ±0.007 


Nemec et al 1988 



Figure 7. The variation of (a) log (Pab) and (b) (A) with [Fe/H], 
taken from the compilation of Djorgovski (1993). Solid symbols 
indicate clusters with 7 or more RRab variable stars. 

An examination of the average periods in Fig. |^ also in- 
dicates that while there does seem to be a linear relation in 
the total sample, if the Oosterhoff groups are considered sep- 
arately the slope is significantly smaller. We have included 
only those clusters having 7 or more RRab stars. We find 

logP= (-0.019 ±0.010)[Fe/H] + (-0.281 ± 0.015) 

for 20 Oosterhoff I clusters, and 

logP = (-0.028 ±0.015)[Fe/H] + (-0.251 ± 0.031) 

for 13 Oosterhoff II clusters. The slopes are consistent with 
each other, and somewhat shallower than predicted by evo- 
lutionary models (—0.05; Lee, Demarque, & Zinn 1990). For 
the two groups, we find average periods of 0.556 and 0.643 
days, respectively. The ofTset in periods between the two 
Oosterhoff groups almost exactly corresponds to the dif- 



ference between the average A values, and so is probably 
not due to differences in mean temperature of the variables. 
From van Albada & Baker's relation, this indicates that 
there is likely to be a difference in the mean mass and/or 
luminosity of the RR Lyraes in the two groups. 

While the Oo II clusters with the bluest HBs show 
large scatter in (A) as expected for small samples of RR 
Lyraes, the clusters with redder HBs (and generally, with 
larger numbers of HB stars) - most notably M68, but also 
M15, M53, NGC 5053, NGC 5466, and possibly NGC 2419 - 
have the same value to within the errors. What also makes 
these clusters unusual is that if they indeed had high he- 
lium abundances as their higher (A) values might indicate 
(whether primordial or due to the action of a deep mixing 
mechanism), all of their HB star distributions would be ex- 
pected to be bluer. For these clusters we are left with the 
possibilities that either high helium is acting against an even 
stronger third parameter, or that helium does not vary and 
there is a different second parameter. 

In short, it would be wise to discard most of the Oo II 
clusters since they are prone to small number statistics, and 
so are overly sensitive to the HB type, making them useless 
as helium abundance tracers. However, the variable-rich Oo 
II clusters stand out in the constancy of their (^4) values as a 
group. Based on these data alone a high helium abundance 
could be a possible explanation. However, the data from 
other indicators show no evidence to support this. It is clear 
that some other factor is necessary to completely explain 
the large numbers (a total of 104 known for the 6 clusters) 
of RR Lyraes in these clusters, but the question is beyond 
the scope of this paper. One result of this discussion is that 
there appears to be no relationship between the Oosterhoff 
dichotomy and the second parameter problem, given the lack 
of variation in (A) with HB type. 
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In the range —1.8 < [Fe/H] < —1.55, there are clusters 
of both Oosterhoff types with at least moderate numbers of 
variables. Three of the Oo I clusters (M3, NGC 3201, and 
NGC 7006) are among the most populated with RR Lyraes, 
and have moderately blue HB types. The Oo II clusters (M2, 
M9, M22, and NGC 5286) are among the bluest HB types 
in the group. The four Magellanic Cloud clusters that have 
{A) between the two Oosterhoff groups also have [Fe/H] ~ 
— 1.8 and have HB types falling between the majority of the 
clusters in the the two groups. 

The difference of 0.06 in A between Oosterhoff I and II 
clusters would require a difference of 0.05 in initial helium 
abundance. A sudden change of that magnitude is unlikely 
to have occurred in Galactic chemical evolution. Enrichment 
of the envelopes of the RR Lyraes by deep mixing on the 
upper RGB could potentially be invoked to explain the A 
measurements. In order to explain observed variations in 
the abundances of species of aluminum and magnesium in 
several clusters, surface material must have mixed through 
regions where hydrogen has at least partly burned to helium 
in the shell source. Sweigart (1997) finds that an increase in 
envelope helium abundance by 0.04 could explain the pe- 
riod shift difference between MS and M15. (Period shifts 
go like 0.84 times the difference in A.) If so, the Oosterhoff 
dichotomy would imply that the driving mechanism only 
existed in one group or the other (in the case of stellar ro- 
tation driving circulation, Oo II clusters would be required 
to have higher average rotation). However, significant vari- 
ations in oxygen abundance are seen in giants in clusters of 
both Oo I (M3, Kraft et al. 1992; M5, Sneden et al. 1992) 
and Oo II (M15, Sneden et al. 1997; M92, Sneden et al. 1991) 
groups, making this explanation unlikely. Among the Oo I 
clusters, the A indicator shows no evidence of significant 
helium abundance variations. 

We have also examined old Magellanic Cloud globular 
clusters and Local Group dwarf spheroidals (dSph) having 
RR Lyrae observations in the literature. Several of the dSph 
galaxies have populations that are old enough to make a 
comparison with globular clusters useful, though composi- 
tion variations have been seen in Draco (Lehnert et al. 1992; 
Shetrone, Bolte, & Stetson 1998) and Sextans (Suntzeff et al. 
1993). In spite of the composition variations, Draco's vari- 
ables have A values that are strongly peaked near its {A) , 
giving additional evidence for lack of a significant metallicity 
dependence in A. 

Ursa Minor, the dSph with the largest mean period of 
those we have data for, falls unequivocally among Oo II clus- 
ters. Carina, Draco, and Sextans have mean periods that put 
them just on the Oo II side of the period gap between the 
Oosterhoff groups. However, only Carina and Ursa Minor 
have (A) values consistent with Oo II clusters. (Our esti- 
mated Rhb value for Carina is uncertain to probably ±0.2, 
but this does not impact the analysis. The value was esti- 
mated from the appearance of the old population HB in the 
CMD of Smecker Hane et al. 1994.) Draco and Sextans have 
{A) values slightly higher than the average for Oo I clus- 
ters. Draco is unusual due to being very metal poor with 
an unusually red morphology. To check if this affected the 
derived {A) values, we derived the following relation from 
the sample of 39 globular clusters (Milky Way and LMC) 
having more than 7 RR Lyrae stars; 



Figure 8. Residuals for {logTun) (in the sense of calculated mi- 
nus observed) as a function of [Fe/H] and Rhb for the globular 
clusters and dwarf spheroidals having more than 7 RR Lyraes 
with good data. 

(logTe//) = (3.8394 ± 0.0012) -I- (0.0100 ± 0.0008) [Fe/H] 

-f(0.0032 ±0.0005)_RHs. 

The rms residual for the globular clusters was 0.0025. The 
dwarf spheroidals do not show good evidence for deviat- 
ing from this rough relation, as shown in Fig. ^. Thus, the 
metallicity and HB morphology seem to be able to account 
for the low temperatures of the RR Lyraes, with no obvi- 
ous distiinction between Oosterhoff groups. As was found 
for the reddest Galactic Oo II clusters, helium cannot be 
the main factor determining the HB morphology since in- 
creased envelope helium would drive the morphology much 
bluer. 

The slope derived from the 5 RR Lyrae- rich Magellanic 
cloud clusters is 

dA/d[Fe/E] = -0.141 ± 0.019. 

This is significantly higher than the slope of the (Oo I + Oo 
II) Galactic samples, but these clusters might also be more 
profitably put into separate Oosterhoff classes. According to 
the mean periods, the SMC cluster NGC 121 and the LMC 
cluster Reticulum are Oo I clusters, and NGC 1786 and 1841 
are Oo II. (We have estimated Rhb = 0.9 ± 0.1 for NGC 
1786 from the statistically subtracted CMD of Brocato et al. 
1996c.) The remaining four LMC clusters (NGC 1466, 1835, 
2210, and 2257) fall between the two groups in Figs. ^, while 
having mean periods at the high end of the Oo I group. 
It may be suggestive that all four of these clusters have 
[Fe/H]fti —1.8, which appears to be at the metal-poor end 
of the metallicity distribution for Oo I clusters in the Milky 
Way. A high helium abundance could be the cause only if 
helium enrichment occurs by a process that did not enrich 
the heavy elements in these clusters. 



3 COMPARISONS 

As stated previously, we have examined these three helium 
indicators in an attempt to see if we could determine with 
greater confidence whether any of the Galactic globular clus- 
ters have abnormal helium abundances. It is also clear that 
we are unlikely to derive useful absolute helium abundances. 
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Table 4. Helium abundance comparisons 
ID (NGC/IC) Y{R) 5y(A) 



5Y(A) 



Figure 9. Correlations between helium abundances derived from 
indicators R and a) A, or b) (A) . A indicates LMC clusters. 



We will now see if we can find anomalous relative abun- 
dances under the assumption that abnormal indicator val- 
ues reflect abnormal helium abundances (which may not be 
the case, given the factors that can affect each indicator). 

For R, helium abundances can be computed from equa- 
tion 11 of Buzzoni et al. (1983). For A, we compute a rel- 
ative heliu m a bundance value (5F(A) using the theoretical 
and the derivative dA/dY 



values in § 2.2 and the derivative dA/dY = 5.8 mag (CCC). 
Here, the uncertainties in the absolute metallicity scale pre- 
vent the computation of good absolute helium abundances. 
We have chosen to use CG abundances from high-dispersion 
spectroscopy, supplementing with abundances from RSD for 
NGC 1261, 1851, and 6218 (M12). For (A) , we compute 
relative helium abundances 5Y{A) by computing the differ- 
ence from the average (A) value for the corresponding Oost- 
erhoff group, and then applying the derivative dA/dY — 
1.4 (Sandage 1990b; Bencivenni et al. 1991). Because the 
(A) values do not show significant dependences on metallic- 
ity, we believe this is the best choice. We have chosen not 
to compute 5Y for clusters having fewer than 7 RRab stars 
with B amplitudes. For the Oo I clusters, the average value 
is 1.8629 ± 0.0017, and 1.9271 ± 0.0028 for Oo II clusters. 
The helium values are given in Table ^, and the comparison 
plots are shown in Fig. 

There isn't an obvious correlation in the comparison 
of helium abundances from the R and A or (A) indicators. 
Comparison between abundances from A and (A) is ham- 
pered by the small overlap between the two samples. We 
can alternately look for clusters whose values may indicate 
anomalous helium abundances. 

In the following paragraphs, we briefiy discuss clusters 
with unusual values for some or all of the helium indicators. 

NGC 1851: The value derived from A is below aver- 
age, and the value from R is about average. From previous 
studies of its RR Lyraes (Wehlau et al. 1978, Wehlau et al. 
1982), NGC 1851 appeared to have an {A) closer to those of 
Oosterhoff II clusters. In fact, both the mean period of the 
variables and the ratio of RR Lyrae numbers Nc/Nab axe 
high for an Oosterhoff I cluster, and the bimodality of the 
HB (with a sparsely populated instability strip) has been 
difficult to model using canonical HB models (see Catelan 
et al. 1998b). However, Walker's (1998) re-examination of 
the RR Lyraes using CCD photometry indicates that the 
_B-amplitudes in the photographic studies were systemati- 
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256 


5272 


(M3) 





204 


5286 








5897 







227 


5986 








6093 


(M80) 





172 


6205 


(M13) 





180 


6333 


(M9) 






6656 


(M22) 






Rup 


106 





233 



Milky Way Clusters 

,Trnm 

-0.023 
,-1-0.028 
-0.034 
1-1-0.019 
■0.022 
-1-0.021 
0.025 



-1-0.013 
-0.015 



1-0.027 
-0.032 
-1-0.021 
-0.024 
-1-0.023 
-0.027 



-0.027 
-1-0.020 
-0.023 
-)-0.020 
-0.022 
-1-0.021 
-0.024 
-1-0.016 
-0.018 
-1-0.015 
-0.016 



-1-0.023 
-0.027 
-1-0.010 
-0.011 



-1-0.020 
-0.023 
-1-0.022 
-0.026 



-0.020 
-1-0.028 
-0.034 
-1-0.030 
-0.037 
+0.026 
-0.031 



-0.024 
-1-0.018 
-0.020 
-1-0.026 
-0.031 
-1-0.011 
-0.012 



1-0.029 
-0.035 



-0.020 
-0.023 
-0.027 



0.044 ± 0.021 
0.048 ± 0.012 
-0.046 ± 0.023 



0.047 ±0.016 

0.001 ±0.021 
-0.039 ±0.019 
-0.016 ±0.016 



0.056 ±0.017 
0.027 ±0.010 



0.017 ±0.016 



0.031 ±0.023 



-0.022 ±0.052 
0.002 ±0.010 



-0.007 ±0.009 
0.008 ± 0.009 

0.009 ± 0.009 
-0.001 ±0.009 

-0.002 ±0.012 
-0.009 ± 0.006 

0.005 ± 0.009 
-0.010 ± 0.006 

0.007 ± 0.010 

0.017 ±0.006" 

-0.002 ± 0.005 
0.009 ± 0.005 
-0.001 ±0.004 
-0.002 ±0.012 

0.012 ± 0.004 
-0.032 ±0.008 
0.023 ± 0.014" 

-0.005 ±0.008 

-0.007 ± 0.006 
-0.009 ± 0.006 



-0.002 ±0.003 

0.002 ± 0.003 
-0.009 ±0.011 

0.013 ± 0.014 



0.009 ± 0.007 
-0.004 ±0.010 
0.014 ±0.001 



cally too large. The new data indicates that NGC 1851's 
(A) value is completely consistent with the Oo I average. 
NGC 5927: see NGC 6624. 

M4 (NGC 6121): AU three indicator values for M4 are 
below the averages, although the {A) value is only slightly 
so. The value derived from A can be questioned due to indi- 
cations that there is differential reddening across the cluster 
(Cudworth & Rees 1990), and that Rv = Ay /'E{B - V) 
differs from the most frequently used value (e.g., Dixon & 
Longmore 1993). Because the magnitudes of the HB and 
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Table 4 — continued 

ID (NGC/IC) Y{R) 
4833 
5694 



5824 
6397 

6809 {M55) 
2419 

5024 {M53) 
4590 {M68) 
7099 {M30) 
7078 {M15) 
5466 

6341 (M92) 
5053 

Reticulum 

1835 

2210 

2257 

1466 

1841 

Carina 
Draco 
Sextans 
Ursa Minor 



220+°-''^^ 



U.ZUU_Q 

n 181 +0-016 



190+°°" 
160+°°2l 

U.1DU_Q Q24 

„ r,c:7+0.018 
U-^O'_0,021 

n 9n«+0-025 
245+°-°^^ 

161+°°^° 
U.1D1_Q 



5Y{A) 



0.030 ±0.019 



0.027 ± 0.009 
0.061 ±0.017 
0.014 ±0.009 

0.005 ± 0.009 



Magellanic Clouds Clusters 
+0.030 



0.150 



-0.037 



0.169"* 



a+0.026 
^-0.031 
204+° °" 
l67+° °23 

U.iD(_Q Q26 
n 1CQ+0.018 

U.±00_Q Q2Q 

Dwarf Spheroidal Galaxies 



5Y{A) 
0.017 ±0.014 

-0.005 ± 0.004 



0.004 ± 0.004 
0.001 ±0.005 
-0.004 ±0.010 

-0.007 ±0.006 
-0.004 ±0.013 



-0.006 ± 0.005 
0.026 ±0.011 
0.012 ±0.014 
0.017 ±0.010 
0.009 ± 0.006 
0.005 ± 0.006 

0.002 ± 0.004 
0.018 ±0.002 
0.022 ±0.006 
-0.004 ± 0.005 



Note: " This cluster's (A) value has probably been affected by its 
HB morphology, and so the 5Y value should be regarded as an 
upper limit. 



MS points were taken from studies of different portions of 
the cluster, we miglit expect a systematic error of over 0.01 
in Y. The reddening quoted in Table ^ is a weighted mean 
of previous studies (Dixon & Longmore 1993), although the 
error of ±0.01 is probably still an underestimate. The values 
from R and (A) should be more reliable. 

NGC 6171 (M107): M107 has very low Y{R) and SY{A) 
values, and an average 5Y{A) value. To completely explain 
the low A value by a reddening error, E(_B — V) would have 
to be underestimated by 0.06, which is possible given that 
E{B — ~ 0.3. Similarly, to explain the low R value, one 
would have to invoke an excessively large metallicity error. 

NGC 6229: Like NGC 1851, NGC 6229 also has a bi- 
modal HB with a sparsely populated instability strip. The 
{A) value derived for this cluster is considerably lower than 
the average. As was previously the case with NGC 1851 (see 
above) , the RR Lyrae data is based on fairly old photometry, 
which may be the source of systematic error in the pulsa- 
tional amplitudes. The R value also indicates a low helium 
abundance for a fair-sized sample (~ 100 each of HB and 
RGB stars). NGC 6229 has a rather extended blue HB mor- 
phology for an outer halo clusters, contrary to what would 
be expected from the (^4) and R values (Borissova et al. 
1997). 

M62 (NGC 6266): Like other post-core-collapse clus- 
ters (see NGC 6752 and M30), M62 has a rather large R 
value, although its {A) value is fairly close to the average. 
M62 has heavy differential reddening, which may affect the 
determination of the HB magnitude. 



M92 (NGC 6341): M92's R value seems to be unusually 
large, although examination of more recent wide-field CCD 
photometry indicates that our value is probably too high 
(Bolte & Roman 1999; in preparation). The A value falls a 
little lower than the average. 

NGC 6624: This cluster has the lowest R value of any 
cluster examined here. Because AVbc is sensitive to [M/H] 
at the metal-rich end, and because NGC 6624 is one of the 
most metal rich clusters in our sample, the metal abundance 
is a natural suspect. NGC 6624 appears to be approximately 
0.2 — 0.4 dex more metal-rich than 47 Tuc according to Ca H 
triplet measurements (RSD; Idiart, Thevenin, & de Freitas 
Pacheco 1997; Armandroff & Zinn 1988), but less than 0.1 
dex more metal-rich according to spectrophotometry (Gregg 
1994). If we optimistically take NGC 6624's metaUicity to 
equal that of 47 Tuc, the R value increases to 0.83 ± 0.10, 
which does not alleviate the problem. It was suggested by 
Richtler et al. (1994) that the total metal abundance [M/H] 
for the cluster is considerably lower than its iron abundance 
([Fe/H]= —0.37; ZW) would indicate, based on comparisons 
between photometry and BV92 models. This could be the 
case if Type I supernovae became the primary source of 
heavy element enrichment. If the iron abundance did not 
actually trace the total metallicity for this cluster (and it 
is in fact more metal-poor), this would mean the the dif- 
ferential bolometric correction we used was too large, which 
would help explain the low R value. The cluster NGC 5927 
also shows an extremely low R value consistent with this 
idea (Ca II triplet measurements indicate that NGC 5927 is 
about 0.4 dex more metal-rich than 47 Tuc). However, M69 
and NGC 6496, which are only 0.1 — 0.2 dex more metal 
poor according to [Fe/H] measures, have R values which are 
closer to the average for our sample. 

NGC 6752: This cluster presents a rather large value 
for R from a large sample of stars. The A value is more 
consistent with the average of the globular cluster system, 
although the HB magnitude is notoriously difficult to deter- 
mine due to a complete lack of stars on the instability strip 
and red HB. 

M30 (NGC 7099): In its high R and A values, M30 
presents evidence for a helium abundance enhancement. The 
cluster's reddening is under some dispute though. Reducing 
the reddening to E(i3 — V) = 0.02 would lower the enhance- 
ment computed from A to 0.033, which is consistent with 
the bulk of the other clusters. The R measurement is more 
certain to be unusually high given the relatively large bright 
star samples. 

Ruprecht 106: One concern in photometric analyses of 
Ruprecht 106 has been the uncertainty in its metallicity. 
Francois et al. (1997) have found [Fe/H]= -1.6 ± 0.25 
from high-resolution spectroscopy. The metallicity uncer- 
tainty plays a small role here though because of the small 
dependences of AVbc (at this low metallicity) and A on 
[Fe/H]. The sample of stars for R is relatively small, so Y{R) 
is consistent with the values for other clusters to within the 
errors. Its RR Lyraes have low A values for the metallicity 
of the cluster though, which may be a result of the apparent 
youth of the cluster (Buonanno et al. 1993), and relatively 
high masses of the variables. 
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Table 5. Helium abundance derived from R as a function of 
N = Nrgb+Nhb- 

Sample Size Y Std. Dev. Afcjusters 

Af> 300 0.199 ±0.010 0.036 12 (2 LMC) 

200 < Af < 300 0.193 ± 0.009 0.031 12 
100 < AT < 200 0.201 ± 0.007 0.036 23 (3 LMC) 



Figure 10. Helium abundances versus [Fe/H] (Djorgovski 1993) 
for the three indicators, a) R: • indicates clusters having greater 
than 200 observed HB stars, o indicates clusters having between 
100 and 200, * indicates clusters with less than 50, and A indi- 
cates LMC clusters, b) A. c) (A) . Points have the same meaning 
as in Figs. H and pi 



4 DISCUSSION 

In Fig. we plot the helium abundances derived from the 
three indicators as a function of metallicity. With the current 
dataset, we find that the three helium indicators R, A, and 
A now yield trends of helium abundance as a function of 
metallicity that are consistent with zero to within the errors 
over the range of [Fe/H] sampled. 

There does not seem to be an obvious trend in the Y{R) 
values as a function of metallicity, although this could be 
masked by the considerable scatter in the points. The mea- 
surement error is primarily the result of Poisson errors. The 
situation will probably improve with the careful examina- 
tion of the HB and RGB populations of clusters with the 
laxgest evolved star populations. A total RGB and HB star 
sample in excess of 1000 stars is necessary to reduce the 
error bars for individual clusters to (^{Y) ~ O.OL 

By separating the clusters into three subsamples accord- 
ing to the total numbers of HB and RGB stars they have, we 
derive the mean values listed in Table ^ The standard devi- 
ations in each case are larger than the average measurement 
error for the clusters in each subsample, indicating that the 
scatter in the measurements may be real. To test this idea, 
we did Monte Carlo simulations using the measurement er- 
rors to set Y{R) for each cluster (with the average value for 
each cluster always centered on the mean of the sample), 
and determined the probability of getting a standard devi- 
ation as large as what is observed. For the 17 clusters with 
Nrgb + Nhb > 250, we find a standard deviation of 0.033 
in Y , and an expected standard deviation of 0.018, which re- 
sults in a probability less than 10~* that the measurement 
errors can explain the standard deviation. For the sample 
with Nrgb + Nhb < 200 (having larger Poisson errors), 
the corresponding probability is 0.06. 

The derived mean abundance is considerably lower than 
the expected primordial value Yp = 0.23. Because we believe 
we have removed the most important (known) systematic 
errors in the R measurements, we should ask what might 
cause the Y{R) values to be low. Brocato, Castellani, & 
Villante (1998) suggest that the uncertainties in the rate 



of the ^^C{a,^y^O reaction lead to an uncertainty in the 
derived helium abundance of about 0.02. Another possibility 
is that there is a systematic effect throwing the AVbc values 
off. 

In lieu of systematic errors in AVbc, we can ask what 
physical processes can slow RGB evolution. Langer, Bolte, & 
Sandquist (1999) have suggested that deep mixing processes 
on the upper RGB may result in lengthened or reduced evo- 
lutionary times for those stars, which could modulate the 
resulting R values. This kind of scenario would mean that 
the helium abundances measured by this method (as well as 
those using HB stars in some way) would be affected since 
some of the helium produced by the hydrogen-burning shell 
would be mixed into the envelope of the stars that undergo 
the process. 

We can also attempt to measure the trend in SY{A) as 
a function of [Fe/H]. The primary difficulty is the clusters at 
[Fe/H]~ —1.3 (essentially on the ZW scale) with abnormally 
low values. There are good reasons for removing M4 and 
M12 from consideration due to differential reddening and 
large reddening uncertainty. Along with M4 (which has a 
relatively low A value also), M12, M107, and NGC 1851 
were the clusters for which we had to use metallicities from 
the RHS study since they were not observed by CG. We find 
from clusters with A measurements (excepting M4, M12, 
M107, and NGC 1851) has a linear trend: 

5Y{A) = (0.016 ± 0.007)[Fe/H] + (0.049 ± 0.011). 

The slope has marginal significance, but systematic error 
in the zeropoint of the metallicity scale affects the slope of 
the best-fit line since the 5Y{A) values of high metallicity 
clusters are modified to a larger degree by metallicity shifts 
than low metallicity clusters. If we reduce the metallicity of 
all clusters by 0.21 dex (equivalent to ignoring a elements), 
we derive a best-fit slope of —0.010 ± 0.007 — a change in 
sign. Thus, the dependence of helium abundance on [Fe/H] 
(in addition to the absolute helium abundance) as derived 
from A will not be certain until the absolute metallicity scale 
is improved. 

The data from the A indicator for the two Ooster- 
hoff groups are individually consistent with constant helium 
abundance as a function of [Fe/H]. It is unlikely that the dif- 
ference in {A) between the two groups is due to a difference 
in helium abundance since a similar jump does not appear 
in the other indicators. 

To explain the apparent difference in A values between 
Oosterhoff groups, there must either be some combination 
of a mean luminosity or a mean mass difference between 
the RR Lyrae stars in the two groups. If the difference were 
entirely due to a difference in mean mass, there would be no 
observable effect in either the A or i? indicators. If there was 
only a difference in the mean luminosity (the Oo II clusters 
RR Lyraes would have to be about 15% more luminous than 
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Oo I variables), Oo II clusters should have R values that are 
13% smaller (since R requires the HB magnitude to compute 
the faint limit of the RGB sample), and A values that are 
higher by about 0.15 mag. Prom 15 Oo I clusters having 
more than 5 variables as well as R values, we find R — 
1.100 ± 0.033. For the 7 Oo II clusters meeting the same 
criteria, we find R = 1.139 ± 0.060. The values for the two 
groups are consistent to within the errors. For A, we find 
5Y{A) = 0.0063 ± 0.0054 and 0.0153 ± 0.005 for 7 Oo I and 
3 Oo II clusters respectively. So again, wc have no evidence 
for helium abundance differences between the groups, but 
the small numbers of clusters with A values makes this a 
weak comparison. 

The difference in mean A value between the groups cor- 
responds to a difi^erence in the mean RR Lyrae mass of ap- 
proximately 20%. This is not so large as to make it unrea- 
sonable that there might be differences in the amount of 
mass loss at the tip of the RGB between the two groups. 
The observations require the Oosterhoff II clusters to have 
variables of lower mass. This goes in a direction opposite 
what is needed to explain the HB morphologies of at least 
some of the clusters. 

Clearly this does not speak to the exact cause of the 
Oosterhoff dichotomy, but it can give a little guidance on 
the details of how the dichotomy is brought about. 



5 CONCLUSIONS 

We begin this section by summarizing what we consider to 
be the most important results of the surveys tabulated here. 

With the corrected differential bolometric corrections 
(larger for more metal-rich populations), we now see no sig- 
nificant evidence for variation in the indicator f? as a func- 
tion of metallicity or horizontal branch type. In particular, 
metal-rich ([Fo/H] < —1) Galactic globular clusters now are 
more consistent with the mean, and Galactic bulge fields 
(Minniti 1995) are also likely to show lower helium abun- 
dances. However, only a handful of clusters have helium 
values Y{R) consistent with the favored primordial value 
Yp ~ 0.23. There is evidence that there is real scatter in 
the R values, over and above what can be chalked up to 
measurement errors. Clusters with anomalously high (M30, 
M62, and NGC 6752) and low (M68, M107, NGC 5927, NGC 
6229, and NGC 6624) values appear over a range of metal- 
licity. 

For A, we find that virtually all of the Oosterhoff I clus- 
ters (and the LMC cluster Reticulum) have values consistent 
with a constant helium abundance. The Oosterhoff II clus- 
ters with the largest number of RR Lyrae variables also ap- 
pear to have a constant (A) value (as do the Ursa Minor 
and Carina dwaxf spheroidals, and the LMC globular clus- 
ter NGC 1841). The remaining LMC globular clusters, and 
the dwarf spheroidals Draco and Sextans have (A) values 
slightly above the average for the Oo I group. For the richest 
Oo II systems, evolution is unable to explain the numbers 
of RR Lyrae variables, and a systematic offset in helium 
abundance also seems unlikely given the evidence from the 
other two indicators. The lack of any correlation between 
HB type and (A) is indication that the second parameter 
problem is completely independent from the OosterhoflF di- 
chotomy. The A measurements give us the best constraint 



on the chemical evolution parameter AY/AZ. Using the la 
error bar on the slope of the A - [Fe/H] relation, we get a 
limit |Ay/AZ| ^ 10, which is consistent with measurements 
from extragalactic HII regions (Pagel et al. 1992). 

For A, our results arc consistent with constant helium 
abundance across the range of metallicity sampled. The ab- 
solute value of the helium abundance as well as the exact 
value of dy/d([Fe/H]) depend on the metallicity scale used. 
The most discrepant clusters can potentially be explained 
by errors in the cluster reddenings. 

Comparing data from different indicators, we find that 
the mean trends are all consistent with constant helium 
abundance for metallicities [Fe/H] ^ —0.7. Wc have exam- 
ined the three indicators so as to use the information to ci- 
ther bolster or dispute claims of unusual helium abundance 
from just one indicator. In examining the clusters with he- 
lium abundances from more than one indicator, we have not 
found convincing evidence that any have abnormal helium 
abundances. Systematic effects clearly appear to varying de- 
grees in the data for all three indicators, but we have not 
been able to determine the cause in all cases. 

We must note that none of the three indicators we have 
tabulated has much data covering the most metal-rich clus- 
ters where evidence of helium enrichment may still reside. 
In general, the photometry for these clusters is most subject 
to field star contamination and heavy reddening, making in- 
terpretation difficult. For R, there is the additional problem 
that the red HB begins to overlap the RGB in the CMD, 
making clean values impossible. The redness of the HB also 
tends to preclude the possibility of RR Lyrae stars, and 
hence the possibility of computing A. Reddening and diffi- 
culties in finding the true HB level makes A values unlikely 
without considerable work. Other means must be devised to 
determine good helium abundances for these clusters. 
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